addpath('npy-matlab-master\npy-matlab')
A1 = readNPY('map_LD_633nm_30K_NiI2_2_x.npy')
A1=A1*1e5/2
surf(A1,'edgecolor','none')
set(gca,'xlim',[1 21],'ylim',[1 21])
view(2)
%colormap(winter)
caxis([-4.5,4.5])
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
xticks({})
yticks({})
colorbar('Ticks',[-4:2:4],'TickLabels',{})
set(gcf,'Position',[0, 0, 380, 300])